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Two of recent progress in lattice QCD approach to nuclear force are reported, (i) Tensor force 
from quenched lattice QCD: By truncating the derivative expansion of inter-nucleon potential to 
the strictly local terms, we obtain central force Vc('") and tensor force Vi{r) separately from s- 
wave and d-wave components of Bethe-Salpeter wave function for two nucleon state with J^ = 
1 + . Numerical calculation is performed with quenched QCD on 32^* lattice using the standard pla- 
quette action at j3 =5.7 with the standard Wilson quark action with K = 0.1640,0.1665,0.1678. 
PreUminary results show that the depths of the resulting tensor force amount to 20 to 40 MeV, 
which is enhanced in the light quark mass region, (ii) Nuclear force from 2+1 flavor QCD with 
PACS-CS gauge configuration: Preliminary full QCD results are obtained by using 2+1 flavor 
gauge configurations generated by PACS-CS collaboration. The resulting potential has the mid- 
ium range attraction of about 30 MeV similar to the preceding quenched calculations. However, 
the repulsive core at short distance is significantly stronger than the corresponding quenched QCD 
result. 
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1. Introduction 

Proton and neutron are the fundamental constituents of atomic nucleus. The interaction among 
them are called as the nuclear force. Together with the structures of nuclei, the nuclear force itself 
has been actively investigated in nuclear physics since its discovery about 75 yeas ago. Phenomeno- 
logically, the properties of the nuclear force are characterized by the attraction at medium distance 
with the OPEP tail [|l|, ^ and the strong repulsive core at short distance. The former is respon- 
sible for nuclei to be bound. The latter is important for various phenomena, such as the stability of 
atomic nuclei, super nova explosions of type II, and the maximum mass of neutron star. While the 
medium to long distance properties are accessible with the meson exchange picture, it is difficult 
to approach the short distance properties. In particular, the origin of the repulsive core has not yet 
been answered beyond phenomenological models. Since the nucleons overlaps at short distance, 
the repulsive core is expected to reflect the internal structure of nucleon. Therefore, QCD is con- 
sidered to be the best tool to reveal the short distance properties of the nuclear force. Indeed, one 
attempted to study the repulsive core with lattice QCD by extending the method of static quark 
potential [Q]. However, the repulsive core was not reproduced from this pioneering work. 

Recently, we have developed a new method to extract the nuclear force between nucleons 
composed of non-static quarks and have found that essential features of the nuclear force are repro- 
duced [^]. In this method, lattice QCD is used to generate the Bethe-Salpeter (BS) wave function 
for a two nucleon state in the center of mass frame: 

WaH {x-y)= ^lim (0 I T [p^ (x, 0«)3 (j, 0)] | A^A^> , (1-1) 

where pa = £abc {u^Cysdh) Uc;a and np = Eabc {uaCy^db) dc-p denote interpolating fields for proton 
and neutron, respectively. Note that this represents an amplitude to find three quarks at x and other 
three quarks at y. At large separation, i.e., \x — y\ —>■ large, it shows a desirable asymptotic behavior, 
which is characterized by the scattering phase shift 5 {k) as 

sin(A:r+5(A;)) 

w(r)-^ — ^ + fors-wave. (1.2) 

kr 

Here, k corresponds to the "asymptotic momentum" measured beyond the range of the interaction. 



which is related to the total energy Pq of the state in Eq. ( |lT| ) as Pq = 2 a /m^ + k'^. The amplitude 



Eq. ( |1.1| ) satisfies the effective Schrodinger equation as 

{y^ + k^)yif{r)=m^ j jVV;vw(r,F') VA(r'). (1.3) 

(For derivation, see Ref. In the r.h.s., Vnn plays the role of the interaction kernel. It is most 
generally non-local, and can be defined to be independent of the total energy of the state [^. After 
the constraints from various symmetris are imposed, the derivative expansion leads us to 

Vnn (r, r) = [Vc (r) + Vt {r)Sn + Vls (r)L . 5 + C»(v2) } 5(r - r') . (1 .4) 

Here, Sn = 3(ai • r)(a2 • r)/r^ - di • 02, L = -ir x V, and S = {di + a2)/2. Vc(r), Vjir) and 
Vls ('■) are referred to as "central force", "tensor force", and "LS force". (Iso-spin dependence of 
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these potentials are understood to be implicit.) These three forces play the most important role in 
the conventional nuclear physics. By truncating Eq. ( |l.4t ) up to the first term, the central potential 
Vc(?";^5'o) and the effective central potential V^{r^S\) have been calculated from quenched lat- 
tice QCD [||]. The resulting potentials possess the repulsive core at short distance as well as the 
attraction at medium distance [^]. Note that, owing to Eq. ([Q]), the method can be extended to be 
more faithful to NN scattering experiments [^. 

The contents are organized as follows. In Sect. |2[, we extend our method to tensor force, which 
is obtained from Schrodinger equation in coupled ^S — Dy partial waves in = 1+. In Sect. |3|, we 
give 2+1 flavor lattice QCD results of nuclear force by using PACS-CS gauge configurations. 



2. Tensor force from quenched lattice QCD 

Tensor force plays an important role in nuclear physics. Together with the repulsive core, it 
has important influences on the structures and the stabilities of nuclei. However, phenomenolog- 
ical determination of tensor force is known to be afflicted with an uncertainty especially at short 
distance due to the existence of centrifugal barrier. 

To obtain tensor force in lattice QCD, we consider Schrodinger equation for = 1+. In this 
case, the wave function has two components, i.e., s-wave component and d-wave component. The 
central force Vc{r) acts separately within these two components. The tensor force Vt(?') provides 
a coupling between these two, and the action of the LS force Vls(?') is restricted within the d-wave 
component. If we keep only the first term in Eq. ([Ol), we can obtain only the effective central force 
V(S*(r;^5'i). By keeping the one more term in Eq. (|r!4|), it is possible to obtain the central force 
Vc (?■)(= Vc{r'^ S\)) and the tensor force VT(r) separately. Note that these two terms give strictly 
local contributions in Eq. ([T^. The effective Schrodinger equation Eq. (13) becomes 



{Ho + Vc{r)+VT{r)Sn)w{r)=EW), (2-1) 

where Hq = — V^/mN, and E = /m^ denotes the non-relativistic energy. For definiteness, we 
restrict ourselves to BS wave function i//^(r) for an state with the azimuthal quantum number M = 0, 
i.e., 

Wali{x-y)= lim {0\T[pa{x,t)n^{y,0)]\NN{f =\+;M = 0)) . (2.2) 
We define projection operators ^ and £2 onto s-wave and d-wave components, respectively, as 

=!^V^«/5(^) = 4l^'^«/^(^"'^)' ^^l-^- (2.3) 

geo 



where O denotes the cubic group, which consists of 24 elements. We multiply ^ and =S to Eq. (2. 1) 



from the left. Since Hq, Vc(r) and VT(r) commute with =^ and £2, Eq. (2.1 ) splits into the following 
two equations as 

HQ[^\if]{f) + Vc{r)[S'wW) + VT{r)[^SnWW) = E[^Y]{r) (2.4) 
Ho[J2Ym + Vc{r)[J2Y]{7) + VT{r)[J2Si2¥]{7) = £[^V^](r). (2.5) 

Note that each of these equations has two Dirac indices a and p. For definiteness, we pick up 
{oc,p) = (1,0) components of these two equations, and solve them for Vc{r) and Vj{r). We arrive 
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at 

Vc(r) =£ + ^([^5i2V^](F)//o[<!^'/^](r)-[^5i2i/^](F)//o[^V^](r)) (2.6) 

^T(r) = ^(-[^V^](r)//o[^V^](r) + [^V^](r)//o[^V^](^)), 

with A(F) = [3^\if\{r)[^Sn^f\ir) - [^\if\{r)[^SnW\^)- Note that, if the d-wave component in 
the wave function vanishes, the first line in Eq. ( p^ ) reduces to Vc{r) = {E — Ho)Y{r) / Yi^^)- 

Numerical calculation is performed with quenched QCD by using Blue Gene/L at KEK. 
The quenched gauge configurations are generated by employing the standard plaquette gauge ac- 
tion at j8 = 5.7, which leads to the lattice spacing a^^ = 1.44 GeV [^. Propagators of quarks 
are generated by employing the standard Wilson quark action with the hopping parameters K = 
0.1640,0.1665,0.1678, which correspond to nty, ~ 731,529,380 MeV, respectively. These calcu- 
lations are performed by using A'conf = 1000,2000,2021 gauge configurations. BS wave functions 
are picked up from the time-slice t — to = 9,S,6, respectively, where the ground state saturation is 
achieved within the error bars. While the periodic boundary condition is imposed along the spatial 
directions, Dirichlet boundary condition is imposed along the temporal direction on the time-slice 
t = 0. Wall source is used on the time-slice t = to = 5 after imposing Coulomb gauge. 

The l.h.s. of Fig. |I| shows BS wave functions of a 7^ = 1+,M = state for 529 MeV. To 
reduce the calculational cost, calculation is restricted to the points on the coordinate axes and their 
nearest neighbors for r > 0.7 fm, whereas all points are calculated for r < 0.7 fm. Note that ^Di 
part of the wave function is multivalued as a function of r, which is due to the angular dependence. 
The r.h.s. of Fig. [l] shows the resulting central force Vc(r) and tensor force VT(r) together with 
effective central force V^^{r) = {E — //o)[^i/A](r)/[^i/A](r). From viewpoint of the meson ex- 
change picture, this shape of VT(r) is expected from a cancellation between the contributions from 
pion and rho meson. Note that V(S^(r) contains the effect of Vt(?') through the 2nd order perturba- 
tion, where one expects that V(^*(r) acquires sufficient attraction to generate the bound deuteron. 
However, we see that the difference between Vc(r) and V^**(r) is quite small, which may be due 
to an artifact of heavy quark mass. Fig. ^ shows the quark mass dependence of tensor force Vj{r). 
We see that tensor force is enhanced in the light quark mass region, which suggests the importance 
of direct lattice QCD calculation employing light quark mass. 

A technical comment is in order. Recall that we used the spin (1,0) component of Eq. (2.6), 
where the 2nd equation vanishes at roc (±l,ibl,±l). This is because the spin (1,0) component 
of d-wave part in the wave function for 7^ = 1+,M = is proportional to the spherical harmonics 
^mZo ( ^ 1 ) °^ 3 cos^ 6 — 1, which vanishes atroc(±l,±l,±l). Although these points are removed 
from the plots, statistical error is accumulated in the neighborhood of these points. (For instance, 
see the points at r ~ 0.5 fm in Figs. |l| and^) It is desirable to improve this in the near future. 



3. Nuclear force from 2+1 flavor lattice QCD 

To compare our results with empirical data, a key role is played by a full QCD calculation 
on a large volume employing a smaller quark mass. PACS-CS collaboration is generating 2-1-1 
flavor gauge configurations on a large volume in significantly light quark mass region [Q. PACS- 
CS gauge configurations are generated by employing Iwasaki gauge action at j3 = 1.90 on 32^ x 
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Figure 1: (o:,j3) = (1,0) part of the s-wave and the d-wave components of BS wave function for a state 
with 7^ = 1 + ,M = (left) and the reconstructed inter-nucleon potentials Vc{r) and V-Y{r) (right). 
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Figure 2: Quark mass dependence of tensor force. 



64 lattice and 0(a)-improved Wilson quark (clover) action with a non-perturbatively improved 
coefficient csw = 1-V15 The lattice scale is determined from uij^, mk and ma inputs leading to 
a^^ = 2.176(31) GeV (a ~ 0.091 fm). Hence, the spatial extension amounts to L = 32a ~ 2.90 fm. 
To calculate nuclear force, we use two series of PACS-CS gauge configurations with (jCud, K"s) = 
(0.13700,0.13640) and (0.13770,0.13640), which correspond to m^^ =^702,296 MeV, respectively. 

To calculate BS wave function, we impose the periodic boundary condition along the spatial 
direction. On the other hand, along the temporal direction, we impose the Dirichlet boundary 
condition on the time-slice t = 32{= Nt/2). We locate the wall source on the time-slice t = with 
Coulomb gauge. Note that the setup is strictly symmetric around the hyper-plane t = aiming at 
doubling the number of data by using the charge conjugation and the time-reversal. (See below.) 

We consider the Euclidean four point correlator of nucleon fields with wall source as 

OV (3.1) 



Gap-a'li'{^,y,t) = {0 T \pa{x,t)np{y,t)p'a,n'p 

where p'^ = L.v,vi£«foc {ua{x)CY5dl (y)) Uc-a{z) and n'^ = Y.x,y;z^obc {uaix)CYsdl (y)) da.^{z) denote 
wall sources for proton and neutron, respectively. Statistical noises are reduced by utilizing the 
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following four symmetries, (i) The spatial translation: Gap-a'p'i^^yj) = Gap-a'^'i^+^^y + ^^t), 
where A denotes an arbitrary 3 dimensional vector, (ii) The cubic group: Ga/3;a'j3'(-^) = Saa{g) 
^PP^S) G^p.^,p,{g^'^x,g-^y,t) Sa'a'ig^'^) Sp,p,{g^^), where g denotes an arbitrary element of the 
cubic group. S{g) denotes the (double-valued) representation matrix of S0(3) in the Dirac bispinor 
space, i.e., S{g) = exp [^OijCOij) with Oij = —^[Yi,Yj] for g = e® with co G so(3). (iii) The spatial re- 
flection: Gap-a'P'{x,y,t) = (Tt))aa iyo)pp G(^p.(^rp,{-x,-y,t) {Yo)a'a' {"fi)) pi p> ■ i^"^) The charge con- 
jugation and time-reversal: Gap-a'p'{^,y,t) = (-Cyo)^^ {-CYo)pp '^Ip-a'pX^^y'-^) i-CYo)a'a' 
{—CY)) pipi- Note that QCD Lagrangian has the charge conjugation symmetry: '^q^^^ = Ccf , 
'loq'ia^^ = q^C, '^Apf€^^ = — A^, where the matrix notation of the gluon field is adopted as 
= A^^T^ for color SU(3) generator This implies the charge conjugation of the compos- 
ite nucleon fields 'loN^^^ = —CN^ and 'loN'ia^^ = —N^C, where N represents an interpolating 
field for proton or neutron. Note that, for correlators in imaginary time, the charge conjugation is 
combined with the time-reversal through the complex conjugation. 

The BS wave function ^f^p ir) for the ground state is obtained from the Euclidean four point 
correlator in the large t region after multiplied by Po,/^/(7^,M) as 



GaP:a'p'{^,y,t)Va'pi{J^,M) = AoYap{x -fj'' ,M)e + • • 



(3.2) 



where V^ipi{J^ ,M) = {G2)a'p' {(^2<yM)a'P' for -^^ = 0+ and 1+, respectively. "• • •" represents 
contributions from excited states, which are exponentially suppressed in the large t region. 
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Figure 3: Full QCD results of (effective) central potentials. The l.h.s shows the results for (fCd^jfCj) = 
(0.13700,0.13640), and the r.h.s. shows the results for (k-„j, k,) = (0.13770,0.13640), where the inset is 
suppressed because of the huge error bar. 

Fig. |(left) shows the full QCD results of the central force Vc(r) = {E - Ho)\lf{x;^ 5o)/i//-(x;' Sq) 
for ^So channel and the effective central force Vf{r,^Si) = {E - Hq)^{x\^ Si) /^{x;'^ Si) for 
channel for {K^d, K,) = (0.13700,0.13640). Here, v/^(f;^ Si ) denotes i//-(x;^Si) = [^V]i^) for nota- 
tional simplicity. Vc(r;' So) and V(S**(r;^ Si) are obtained from BS wave functions on the time-slices 
t = S and t = 9, respectively, where the ground state saturations are achieved within error bars. 
Similar to the quenched results[Q, |8]], they have attractive pockets of about 30 MeV in the medium 
distance, i.e., 0.5 ^ r ^ 1.0 fm. In contrast, the repulsive cores are considerably strong. They 
are by about 10 times as strong as the quenched result for comparable pion mass [^]. There seem 
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to be several possible reasons, (i) A dynamical quark effect, (ii) The action adopted in quenched 
calculation may not be close to the continuum limit. 

Fig. fright) shows the full QCD results of the (effective) central potential for {Kud^Ks) = 
(0.13770,0.13640). These results are obtained from the BS wave functions on the time-slice t = 6, 
where the ground state saturations are achieved within the statistical errors. We do not show the in- 
set because of huge statistical errors. Although it is necessary to improve the statistics significantly 
to reduce the huge error bars, we see that the repulsive cores are again considerably strong. 



4. Summary 

We have presented preliminary lattice QCD results for the tensor force by using quenched 
QCD. We have seen that the tensor force has a large quark mass dependence and is enhanced as the 
quark mass decreases. We have presented preliminary results of the (effective) central force from 
2+\ flavor lattice QCD by using PACS-CS gauge configurations. A remarkable difference from the 
quenched results was found in the strength of the repulsive core. It is interesting to investigate the 
reason, since it may provide us with a key to the origin of the repulsive core, which is one of the 
most important open problems in the nuclear physics. 
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